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Abstract 

Determining if a point is in a polygon or not is used by a lot of ap- 
plications in computer graphics, computer games and geoinformatics. 
Implementing this check is error-prone since there are many special 
cases to be considered. In this paper we present a simple even-odd 
algorithm to solve this problem for complex polygons in linear time. 

1 Introduction 

At a first glance, the point-in-polygon problem seems to be a rather simple 
problem of geometry: given an arbitrary point Q and a closed polygon P the 
question is whether the point lies inside or outside the polygon. There exist 
different algorithms to solve this problem, such as a cell-based algorithm 
|ZK01| . the winding number algorithm [HA01J or the even-odd algorithm 
that is used in this paper. The problem is not as trivial as it seems to be if 
the edges of the polygon can intersect other edges as seen in figure [TJ This 
kind of polygon is also often called a complex ploygon since it can contain 
"holes". 

A lot of special cases have to be considered (e.g. if the point lies on 
an edge) and most of the existing even-odd algorithms either fail one or 
more of these special cases or have to implement some kind of workaround 
|Sch08| . We present an even-odd algorithm that works with all points and 
all polygons with no special cases to be considered. 

2 The even-odd algorithm 

The algorithm itself is fairly simple: cast an infinite ray from the point in 
question and count how many edges of the polygon the ray intersects. 



Figure 1: A self-intersecting polygon 

Definition Let P be a polygon with n vertices Po, Pi, P n -i, P n whereas 
the vertices are sorted in such a way that there exists an edge between Pi 
and Pi+i for every index < i < n and between Po an d P n - The algorithm 
computes if the arbitrary point Q lies either inside or outside of P. 

The steps of the algorithm To ease the computations the algorithm 
defines its own coordinate system with (0|0) = Q by moving the polygon. 
The algorithm then uses the positive x-axis as ray to calculate intersections 
with the edges of P. So the start vertex of the ray is (0|0) and the end 
vertex is (x max \0) where x max is an x value greater than that of any of the 
vertices of P. 

1. first it is determined if Q is equal to any of the vertices of P or lies on 
any of the edges connecting the vertices. If so the result is inside. 

2. a vertex P s that does not lie on the x-axis is searched in the set of 
vertices of P. If no such vertex can be found the result is outside. 

3. beginning from that vertex P s the following steps are repeated until 
all vertices of P have been visited: 

(a) the index s is increased to s + i until the next vertex P s +i not 
lying on the x-axis is found. If the index s + i > n then % is set 
to —s and the search is continued. 

(b) Depending on the course of step (a) one of the following steps is 
taken: 

i. no vertex has been skipped: the line segment from P s to P s+ i 
is intersected with the positive x-axis 

ii. at least one vertex with a positive x-value has been skipped: 
the line segment from P s to P s +i is intersected with the com- 
plete x-axis 
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iii. at least one vertex with a negative x-value has been skipped: 
nothing is done 

(c) P s +i is the starting vertex for the next iteration. 

4. if the count of intersections with the x-axis is even then the result is 
outside, if it is odd the result is inside. 

Proof of correctness To proof that the oven-odd algorithm in general 
is correct one can generalize it to the Jordan Curve Theorem and proof its 
correctness |Tve80| . What has to be proven is that the given algorithm is in 
fact a correct even-odd algorithm, meaning that the count of edges is correct 
under all circumstances. 

A lot of problems emerge when trying to intersect two line segments and 
one of the segments has one or two of its vertices lying on the other segment. 
If the count of intersections with the x-axis are to be counted correctly, then 
each edge of P must either clearly intersect the x-axis or not intersect it at 
all. There must be no case where the starting or end vertex of a line segment 
lies on the line segment it should be intersected with. 

The x-axis is the first line segment to look at, since it is part of all the 
intersections. The start vertex of the x-axis, namely Q, is guaranteed not 
to lie on any edge or to be equal to any vertex of P. If this was the case, 
then the first step of the algorithm would have already returned the correct 
result. The end vertex is guaranteed to have an x value greater than any of 
the vertices of P, so no vertex of P can be equal to it and no edge of P can 
contain it. 

Of course there still exists the problem that one or more vertices of P 
are lying on the x-axis as seen in figure [2] 




Pi 



Figure 2: One of the edges of P lies on the x-axis 

The algorithm deals with this kind of problem by ignoring the vertices 
lying on the x-axis and stepping over them when trying to find an edge to 
intersect. It then creates a new auxiliary edge that it can intersect safely. 
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So starting for example at vertex P\ it would ignore P% and P3 and then 
create a new edge between Pi and P4 as shown in figure [3j 




Figure 3: A new auxiliary edge has been created 

At this point it is clear that none of the edges used to calculate the count 
of intersections with the x-axis has a vertex on the x-axis and is either clearly 
intersecting it or not. What has to be shown is that all created auxiliary 
edges are correct substitutes to calculate the intersections with the x-axis. 

Indeed they are not a correct substitute to intersect the positive x-axis, 
as can be seen in figure [4j Here the auxiliary edge would be the same as the 
edge between P\ and P3 and this edge does clearly not intersect the positive 
x-axis. So the total count of intersected edges of the polygon shown in figure 
[4] would be zero - which is obviously wrong. 

The algorithm actually deals with this problem in step 3.6, by extending 
the ray in that special case where a vertex lying on the positive x-axis has 
been skipped. The new ray is then the complete x-axis and not only the 
positive part of it. As can be seen this would create the desired result in 
the example of figure |4j because the total count of intersected edges would 
then be one. 




Figure 4: The auxiliary edge between Pi and P3 does not intersect the 
positive x-axis 
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The question remains if this method always yields the correct result 
under all possible circumstances. If the skipped vertex lies on the negative 
x-axis then the auxiliary line will not be considered for intersection, since 
none of the original edges could have intersected the positive x-axis. So all 
cases that have to be looked at involve one or more vertices on the positive 
x-axis. The skipped vertices can never change from the positive to the 
negative x-axis since this would have been caught in the first step of the 
algorithm. 

So all auxiliary edges that intersect the x-axis are created from the fol- 
lowing order of vertices: P u which does not lie on the x-axis, P^\ which does 
lie on the positive x-axis and P w which does not lie on the x-axis. Each 
of the vertices P u and P w can lie in one of the four quadrants around the 
point Q and therefore there exist 4*4 = 16 different versions of the auxiliary 
edge that have to be considered. This number can be further reduced to 10 
because six of these versions are created by switching start and end vertex 
and do not have to be considered as a separate case. 

In table [2] all possibilities of the locations of P u and P w are shown as 
well as the desired intersection count and the actual intersection count of 
the algorithm. It is clear by looking at the table that the algorithm satisfies 
the desired result for each possible scenario. Therefore the auxiliary line is 
indeed a correct substitute for the original edges. 

This leads to the conclusion that the provided algorithm can correctly 
determine if there is an even or an odd number of intersections with the 
polygon. 

□ 

Time complexity All n vertices of the polygon are visited once during the 
translation and the first three steps of the algorithm. All other computations 
like step four can be completed within constant time. 

Therefore the time complexity for this algorithm is 0{n). 

□ 

3 Example 

The algorithm is applied to the sample complex polygon P in figure [5j 

1. Q lies neither on any vertex nor edge. 

2. the start vertex P s is Pi in this example. 

3. all intersection and substitution steps can be found in table [TJ 

1 Of course there could be more than just one vertex on the x-axis but they are all 
skipped and the same auxiliary edge is created no matter how many additional vertices 
are on the x-axis. 
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Pi 

Figure 5: Sample complex polygon P 



4. since the count of intersections with the x-axis is odd the result is 
inside. 
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Table 1: Intersection and substitution steps 



4 Implementation 

The translation of P can be done together with the first step of the algorithm 
by moving the vertices by the negative x- and y- values of Q. The geometric 
helper functions to intersect line segments can be found in [Stu99] which are 
an improved version of [Sed92 . 

5 Conclusion 

The implementation of point-in-polygon algorithms for complex polygons is 
error-prone since there are many special cases to be considered. We pre- 
sented a correct even-odd algorithm that prevents special cases by substi- 
tuting a sequence of edges with an auxiliary edge. 

The correctness of the algorithm has been proven for all possible poly- 
gons. 

The algorithm is also time-efficient since it is 0{n) for polygons with n 
vertices. 
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A Line substitution possibilites 
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Table 2: Possible cases when creating an auxiliary line (blue). The values 
ql to q4 correspond to the quadrants of a cartesian coordinate system (ql : 
x>0Ay>0,y-2:x<0Ay>0, ...). Rows that are just gained by switching 
the vertices P u and P w are shaded gray. 
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